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Abstract 

The Hamiltonian and trial function in the BCS theory are improved to test 
the limit of this theory. The Cooper pairs arise from standing electron waves, 
ready to move with atoms, giving high Tc. The Hamiltonian is derived from 
alternating layers of light and heavy atoms, giving a forbidden zone hosting 
no standing wave pairs. The exchange term may force singlet pairs into this 
zone, leaving triplet pairs outside, giving magnetic excitations. If the Fermi 
energy is crossed only by the Cu02 band, then the forbidden zone and triplet 
pairs will vanish, consistent with experimental evidence. 



PACS: 74.20.Fg - BCS theory and its development 
PACS: 74.80.-g - Spatially inhomogeneous structures 
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Recently there has been considerable interest in the original theory of Bardeen, 
Cooper and Schriefer (BCS) [Q. It was used by Nunner, Schmalian and 
Bennemann to explain the isotopic effect of cuprates [Q. It was also used 
by Nozieres and Pistolesi to study pseudogaps [^. Following this line of 
thinking, we adopt the view that cuprates have a Fermi surface and electron- 
phonon interactions. We also adopt the BCS formalism, which is relatively 
simple, due to neglecting damping and retardation, with negligible effects [Q. 

The BCS theory is based on variational evaluation of the Frohlich Hamil- 
tonian, with the Cooper pairs as the trial function ^. With an im- 
proved Hamiltonian and trial function, we may explain further properties 
of cuprates. In order to improve the Hamiltonian, we use a simple model 
of alternating layers of light and heavy atoms (Fig. 1), because in cuprates 
the Cu02 layers are always sandwiched by heavy atoms. For example, in 
La2Cu04 we have a ratio 96/310 when comparing the atomic wt. of the 
Cu02 layer and the two LaO layers. In YBa2Cu307 we have a ratio 281/386 
when comparing the two Cu02 layers (including the Y atom in between) with 
other atoms. In Tl2Ba2Cu06 the ratio is 96/746, which becomes 256/746 and 
416/746 when the number of Cu02 layers are 2 and 3 (including Ca in be- 
tween). In HgBa2Cu04 the ratios are 96/507, 232/507 and 368/507 when 
the number of Cu02 planes are 1, 2 and 3 (including Ca in between) 0. 

The classic Cooper pairs arise from traveling electron waves, which are 
mobile in all directions. However, the resistance of cuprates in the c axis is 
significantly larger than that in the a-b plane most electrons are mobile 
only in that plane. It is reasonable to assume that the atomic layers act as 
potential wells to retain electrons. Being reflected back and forth in potential 
wells, electrons must be in the form of standing waves. Indeed, according to 
energy band calculation, in cuprates the Fermi surfaces join together in the c 
direction |]^, a familiar sign of electron standing waves in that direction. In 
the case of Bragg scattering only some electrons are reflected, giving narrow 
'necks' to join Fermi surfaces. In cuprates virtually all electrons are reflected 
in the c direction, so that the Fermi surface becomes cylindric. In fact, many 
cuprate theories assume that carriers are somehow bound around the Cu02 
plane ||^, i.e. they are standing waves across that plane. We add two Bloch 
functions together to model electrons in cuprates. The result is a traveling 
wave in the a-b plane but a standing wave in the c direction (Fig. 1). The 
Cooper pairs arise from such waves. These pairs concentrate in the layers of 
light atoms (Cu02 planes) as a natural result of our theory. 

First, we explore the prediction power of the BCS theory. The reduced 
Hamiltonian [|T| involves a series of pair generation and destruction operators 
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with the c-number coefficient 



l^(k,q)=X:' 2;^a;,(q)^f(k,q) 



where k and e are electron vector and energy (measured from the Fermi 
surface), q and ui phonon vector and frequency, and / identifies phonon 
branch (excluding transverse phonons, which do not interact with electrons 
in iV-processes) 0. The matrix element 



Mi(k,q) = qi 



hN 


1/2 


2MLJi{q)_ 


Jn 



^k+q,.(r)5V(r)V^k,.(r)c?r (2) 



measures the strength of electron-phonon interaction. Here ip is the electron 
wave function, a spin t or |, Mthe mass of an atom. A?" the number of atoms 
in unit volume, r the coordinates in real space, flo a volume surrounding the 
atom. To its boundary, and 5V(r) = V(r) — V(ro), V being the potential field. 
We define qi as the l-th component of t/q, U being the 3x3 unitary matrix 
found when solving the classical equation of motion for the atom. Mott and 
Jones found matrix elements when Qq is the Wigner-Seitz cell [ITDI. We find 



equation ^ when V(ro) is constant (this defines f2o in a natural manner). 
The BCS self-consistent equation 



is an integral equation of the Cauchy type: ^(k, q) is singular 

We solve equation ^ through iteration |]1T|. With a proper ffist ap- 
proximation, we may expect reasonable accuracy after just one iteration. 
We use free electron energy to evaluate equation [^, where the denomina- 
tor turns out to be Aepe^lSf — (C + cos 6*)^], ep = (/i^/2m)|kp is the Fermi 
energy (we study |k| near the Fermi surface), eq = (/i^/2m)|qp, 5f = 
{m/2)vf /ep, vi = co';(q)/q| the sound velocity. In the Debye approxima- 
tion 6i = {Z/IGY^^Qd/Qf ~ 10^'^ in all superconducting metals, where Z 
is the valency, and Qp are the Debye and Fermi temperatures. There- 
fore V^(k, q) > (condition to have an energy gap) holds in equation |I] only 
when ( + cos9 ~ 0, C = |q|/|2k|, 9 being the angle between k and q, so 
that |k + qp = |kp + |qp + 2|k||q| cos6' ^ |kp. Thus e(k q) ^ e(k), i.e. 
electrons changes direction but not energy in scattering, which is used as our 
ffist approximation (also used to study metal resistivity) [jlO[- The use of 
free electron energy implies a spherical Fermi surface and hence an isotropic 
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energy gap, so that A(k + q) = A(|k + q|) ^ A(|k|) = A(k). These lead 
through equation ^ to: 

Etrv = E ^(k, q) = jr^Vpnv^ (4) 

where Etrv = [A2(k) + e^ik)^/^ at T = 0, TRV standing for travelhng 
wave, e and n are electron charge and density, ks the Boltzmann constant, p 
the resistivity at temperature Tp, v = ksQu /T^ko the Debye sound velocity, 
ko being the phonon cut-off wavenumber, and 

1 /.(4Z)-l/3 /•2J/- , .{4Z)-i/3 

ri=- ^'(^)r4^// (5) 

TT JO i — C^/ JO 

Here = 3(xcosx — sinx)/x'^ is the overlap integral function, x = 

SMa^/^Z^/^C, a = NQo/^ the fraction of Qq in a primitive cell, and f2 
the unit volume. We assume |q|/| 2k| = C < (4Z)-i/3 < 1^ because equa- 
tion |I| arises from a canonical transformation ^ , where operator commuta- 
tion requires q 7^ ±2k. In first iteration, the pair occupancy varies linearly if 
—Etrv < e(k) < Etrv, otherwise equation |]has improper solutions (occu- 
pancy < or > 1). This justifies the BCS approach to integrate equation | 
only in a thin layer across the Fermi surface This surface does not have 
to be spherical, so long as kn « |k| (integration area small). 

In order to find Tc we follow BCS to minimize the free energy of the 
electron-phonon system. Evaluating the result via iteration, we find 

E = Etrv tanh{E/2kBT) (6) 

where E = [A2(k) + e^k)^/^ at T > 0. Since tanh(E/2A;BT) < E/2kBT, 
we have T < ETRvf^ks and hence 2Etrv /ksTc = 4. It is easy to prove, 
by direct substitution, that E from equation ^ is also the exact solution of 
equation 3.27 in (BCS self-consistent equation for T > 0), provided that 
hujN{0)V = Etrv in that equation. Clearly, E is not a function of k, the 
so-called gap parameter A(k) is, contrary to general perception. This is 
justified: according to BCS it is E that measures the energy gap, A is just 
an approximation, i.e. e ~ near the Fermi surface giving E = (A^-|-e^)^/^ ~ 
A [|l|. What arises from experiment is actually 2Etrv IksTc, whose value 
may deviate from 4 for reasons other than huj/kBTc >> 1 (weak coupling). 
Since E is constant, any Cooper pairs are equally likely to be excited. This 
is also justified: M.i (measuring the strength of electron-phonon interaction) 
in equation |^ varies little across the Fermi surface. 

According to equations ^ and || a good superconductor must have numer- 
ous free electrons (large n) scattered frequently by atoms (large p) moving 
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quickly to facilitate pairing (large v). The factor rj arises when the sum- 
mation over q in equation ^ is replaced by an integration over (47r/3)fc|, 



(volume of the first Brillouin zone), which exists in the sense of the Cauchy 



principal value (used by Kuper to verify the BCS theory) |TT|, 0, i.e. posi- 
tive and negative contributions of ^(k, q), if finite, are cancelled on a series 
of spherical surface, the singular point ignored. We are entitled to do so, 
because equation |l| is defined on a grid of k and q, which may not be in 
precise combinations to let V^(k, q) = oo. We can also avoid such com- 
binations by suppressing a few phonons with little physical consequence. 
This principal value varies little among phonon branches, allowing us to use 
J2i Qi = q*t/*f/q = |qp {U unitary) to find the numerator in equation ^. We 
use the expression for metal resistivity to calibrate 6V, and this leads to the 
denominator in equation |^ [|l^]. When a = 1, equation || yields 2Etrv = 2.2, 
18 and 27 for Cd, Ta and Nb (1.5, 14 and 30.5 experimentally, in 10~^eV), 
which are of the right order, although over and under-estimtions are possible. 
On average equation ^ yields 2Etrv = 15.7 for Zn, Cd, Hg, Al, Ga, Tl, Sn, 
Pb, V, Nb, Ta and Mo (11.3 experimentally). 

Now consider a crystal of alternating layers of light and heavy atoms 
(Fig. 1). For Cooper pairs of traveling electron waves, equations ^ and ^ are 
still valid. The derivation is straightforward in principle but involved tech- 
nically. In equation |^ the upper limit of integration is replaced by (8Z)~^^^: 
the first Brillouin zone is smaller the larger the primitive cell. We also have 
a = 2NQo/Q (assuming Qq invariant in the cell) and Z averaged over dif- 
ferent atoms. In order to model electrons in cuprates, we notice that, in 
principle, electrons of any configuration can be expanded into Fourier series 
in terms of plane waves. These series can be shortened, when the plane waves 
are replaced by waves resembling more closely the actual configuration of the 
electrons in the crystal, as is done by various sophisticated methods to cal- 
culate the electron band structure 0. For simplicity, we consider series of 
superpositions of just two Bloch functions: 

ip^lir) oc exp[i{kr,x + kyy)] cos{'kz/c) (7) 

'(/'^(r) cx exp[z(fc^a; + kyu)] sin(7r2;/c) (8) 

where k = xA;^ + yky is a 2D wavevector. Both il)^^"^ and are traveling 
waves in the a-h plane. On the other hand, if)^^^ is confined in layer 1 (its 
anti-nodes are in that layer), whereas t/'*^^-* is confined in layer 2 (Fig. 1). 
Apparently, electrons overlap with those in the neighboring layer, but not 
with those in the next neighboring layer (there is a node in between, see Fig. 
1). These appear to be reasonable as first order approximation. We assume 
that both ip^^'^ and ip'^'^'^ are close to the Fermi surface. This is true at least 



5 



in some cuprates: the chain band crosses ep in YBa2Cu307, so does the TIO 
band in Tl2Ba2Cu06, the BiO band in Bi2Sr2CaCu208, etc. 0. 

We use standing electron waves to build Cooper pairs. The formalism 
parallels that for the classic Cooper pairs: in second quantization the basis 
states do not have to be plane waves [0, |l3l. This is necessary, because super- 



conductivity is a second order process. Unless the configuration of (single) 
electrons is chosen properly (a first order process), the energy gap cannot be 
maximized. In equation V is replaced by Vij = 1,2), / runs from 1 to 
6, and A4i is replaced by A^f"*^ to link ip'^^^ and Both intra {i = j) and 
inter-layer {i ^ j) couplings are possible, because electrons in neighboring 
layers overlap (Fig. 1) . If we use c/ 2 to replace c in equations ^ and H, then we 
are out of the first Brillouin zone. As a result, t/'*^^-* and tjj^'^^ become identical, 
giving no inter-layer coupling and reduced energy gap, as will be shown. If 
we use 2c, 3c, ... to replace c, then E12 in Appendix becomes smaller. This 
weakens the inter-layer coupling, a choice not preferred. Indeed, the classic 
Cooper pair (spin |, | and wavevector k, — k) is also a choice (i.e. a trial 
function) to maximize the energy gap. 

Letting f^^ and f^^ be the over-all probability of excitations in layer 1 

and2,wehave-4fcBEk[(tfV2 + /f /2) + (l-tfV2-/f /2)ln(l-/«/2- 
/i^^''/2)] as the entropy of the pair ensemble, which is from the consideration 
that pairs in either layer may have the same energy, so that thermodynami- 
cally they fall into the same group of entities, i.e. there can be 4 electrons at 
the same energy level, giving the degeneracy factor 4. Note that in cuprates 
bands of both the light and heavy layers may cross ep so that the above 
degeneracy is possible. Minimizing the free energy of the pair ensemble, we 
find through iteration two equations with the solution 



2/k = Tr2 ^r-^keT In (9) 

1 g ^(2) _ Eu - En 2 - fi^^ - fl^> 

E12 - ^11^22 /k + /k 

Here Eij = Z]q^j(k, q) and the summation over q = xg^. -|- yqy is in 2D, 
i.e. the standing electron waves emit and absorb phonons in 2D. Since Eij = 
N~^J2Eij when the summation is over (not an argument of Eij), 
being the number of q^ in the first Brillouin zone, we integrate E^j in 3D 
(in the sense of the Cauchy principal value) for convenience. Equations § 
and apply to pure intra-layer couplings when £"12 = 0, and pure inter- 
layer couplings when En = and E22 = 0. Both lead to an energy gap, 
because E^j > always holds (Appendix). 
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The inter and intra-layer couplings are in competition. If none dominates, 
then both are suppressed. Since Mi < M2, we have E22 < En (Appendix). 
When E22 < E12 < {EiiEi2Y^'^ , the two sides of equations ^ have opposite 
signs: f^^ = 1/2 and T = must hold to give f^^ = 1/2 via equation 
Similarly T = when (^11^22)^/^ < ^12 < ^11 • Clearly E22 < E^ < Eu 
(inter-layer coupling weaker than the coupling in layer 1 but stronger than 
that in layer 2) is a forbidden zone hosting no standing wave pairs. Outside 
this zone, equations | and |10| can be added together to recover equation 
with E replaced by 1 — f^^ — f^^ and Etrv replaced by 

E "^12 — -^11-^22 (11) 

E12 — {Ell + -E22)/2 

STD standing for standing wave. If all Cooper pairs are excited, then 
fL^\fk ^ ^ 1/2 giving 2EsTD/kBT, = 4 to estimte T,. 

In Fig. 1 Estd/Etrv > 5.57 when Eu < E12. This large ratio arises 
from equation |lT], where Eu Eu leads to Estd — 2£'ii: the energy gap 
is larger the stronger the inter-layer coupling. At this point, equation H yields 

('2') 

/^"^ = 1/2: all the pairs in layer 2 are excited, apparently draining much of 
the excitation energy. Pairs in layer 1 (light atoms) are more or less left alone: 
superconducting carriers are in the CUO2 layers. Furthermore, Eu > Etrv 
holds as a result of the symmetry of the standing waves (with respect to a-b 
planes, rather to sites of atoms) reflecting the fact that bound electrons are 
readier to move with the atoms. In Fig. 1 Eu < E12 for Z < 0.1314. This 
small valency arises, because on average phonons have smaller |q| to pair 
electrons on a smaller Fermi sphere, so that Ai[^^\(X |q|) < M.\^'^\(x 2tx jc) 
holds to give Eu < -^12- Assuming Estd/Etrv = 5.57, we have ~ 130K 
when Etrv = 40.2 x lO'^eV (40 x lO'^eV for NbgGe). In Ba2YCu307 we 
have ko = (Gn^N/ny/^ = 6.99 x lO^m'^, 9^ ^ 400K, v = kB^D/hko ~ 
7.49 X lO^ms^^ p = 70 - 550 x 10"^fim {a-h plane) and n ^ 6 x lO^^m'^ 
from infrared reflectivity [1^. Taking surface values, these lead through 
equation ^ to Etrv = 10 — 77 x lO'^^eV (77 = 1). A point to notice: p 
is for traveling waves in equation ^ Although standing electron waves are 
easier to be scattered, giving larger p, this is more or less compensated by 
the weak coupling at their nodes. It is interesting that, in Ba2YCu307, 
= 5.76 X lO^^rn-^, so that Z = {l/13)n/N ^ 0.08 i). 
We may have spin singlet pairs in the forbidden zone, triplets outside. 
Specifically, while £"12 < Eu holds for singlet pairs, £^12 > £^11 may hold 
for triplet pairs, i.e. Eu declines faster when the pair symmetry changes, 
due to the stronger effect of the exchange term on Eu, which is related 
to intra-layer coupling, where electron waves overlap to a greater extent. 
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Outside the forbidden zone, Estd changes ~10% when Z drops just 0.001 
from 0.1314. Across the zone border, Zstd changes more dramatically (Fig. 
1). In both cases doping may obscure the isotopic effect. The Fermi surface 
of cuprates is not strictly cylindric 0: classic Cooper pairs may arise to give 
superconductivity inside the forbidden zone. Since Etrv changes slowly with 
Z (Fig. 1), the isotopic effect will be more apparent. Indeed, in cuprates the 



isotopic effect is minimum when Tc peaks with proper doping . Although 
travelling electron wave pairs (giving Etrv) cannot compete with standing 
wave pairs (giving Estd > Etrv) in the a-b plane, they can move in the c axis 
(standing waves cannot) to make the Knight shift complicated ]TB| , |T^, and 
give energy gap anisotropy and pair symmetry anisotropy (well documented 



for cuprates) [TS 



In conclusion, rather surprisingly and significantly, the BCS theory may 
play a major role to explain the properties of cuprates. We show that this 
theory can be used to calculate from first principles, the first time to our 
knowledge, giving Tc ~ 130K for cuprates. The impression that Tc is low 
in the BCS theory arises from McMillan's calculation, where the Tc of an 
alloy family is clamped to that of the parent metal, assumed to be a natural 
element of low Tc |]2^. Another worry about the BCS theory is the Migdal 
instability which, according to Waldram, may not set in at 130K |1^, p5| . 
We also show that in a complex system like cuprates the microscopic physics 
may manifest itself as a paradox: the exchange term may drive the singlet 
pairs into the forbidden zone, leaving the triplet pairs outside. Indeed neu- 
tron scattering does exhibit a magnetic peak below Tc from YBa2Cu308 and 
Bi2Sr2CaCu208 [^, and the references therein]. In addition to adding an- 
other explanation to the already long list of explanations, our theory has 
a specific experimental basis. In the above cuprates ei? is crossed by both 
the Cu02 and heavy atom layer bands [0]: our theory applicable. On the 
other hand, magnetic excitations are absent in La2_a;Sr2.Cu04 where 
ep is crossed only by the Cu02 band [0. Therefore we have to assume 
Ml = M2: heavy atom layers are not involved in the electron-phonon in- 
teraction near the Fermi surface. As a result, both the forbidden zone and 
triple pairs vanish. It appears worthwhile to search magnetic excitations from 
e.g. HgBa2Ca2Cu308 and Tl2Ba2Cu06, as example and counter-example of 
cuprates, where ep is crossed only by the CUO2 band |]7|, p8|. 
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Appendix 



For a spherical Fermi surface 



Eu = C [ 
Jo 



(8Z)-V3 ^ Uji^B)'' [A-^Bf] CdC 




(8Z)-V3 ^2 1 0.65 a3 C2^C 

2 [Mi ^ M2J ^273 ^(1 -(2)2 



yl and 5 are values of the overlap integral function F{x), with x = 3.84a^/^Z^/^ 
C(l - Cy/^ and 3.84ai/3[^2/3^2(i _ ^2^ ^ o.65(a/c)^/3]i/2^ respectively, and 

C = 6Zma'^{6Vy /ep- If Mi and M2 are interchanged, then En is turned 
into £'22- For Cooper pairs of traveling waves 5 = and X = 3Ma^^^Z^/^(, 
so that Ell is reduced to Etrv in equation H (C expressed in n. p and f ). 
The factor (1 — C^)^^ in -E'li (or Etrv) and £'12 is related to the Cauchy 
principal value, which starts to fail when = |q|/|2k| — >■ 1, where the canon- 
ical transformation also fails. However, in Fig. 1 Etrv shows little sign of 
divergence when Z > 0.13. Another factor (1 — C^)^^ in -£'12 is from the de- 
nominator (oc |q|2) in equation |I], which is cancelled in En [Aij^ oc |q|) but 
not in E12 (oc 27r/c), so that in Fig. 1 Estd turns upwards when Z 0.13. 

Figure legend 

Fig. 1 Crystal of light (wt. Mi in layer 1, open circles) and heavy atoms (wt. 
M2 in layer 2), a,b and c are lattice constants, the anti-nodes of ip^^^ and ip^^^ 
are in layer 1 and 2, respectively; Estd (solid line, a.u.) and Etrv (broken 
line, a.u.) are found when M1/M2 = 0.7 and a = b = 0.5c; Z < 0.1314 
{Ell < E12) and Z > 0.1361 (£^12 < £22) border the forbidden zone. 
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